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Abstract 

We study the Langevin dynamics of a heteropolymer by means of 
a mode-coupling approximation scheme, giving rise to a set of cou- 
YY ■ pled integro-differential equations relating the response and correla- 

^^ i tion functions. The analysis shows that there is a regime at low tem- 

perature characterized by out-of-equilibrium dynamics, with violation 
of time-translational invariance and of the fluctuation-dissipation the- 
orem. The onset of ageing dynamics at low temperatures gives new 
C^ I insight into the nature of the slow dynamics of a disordered poly- 

mer. We also introduce a renormalization-group treatment of our 
FTj . mode-coupling equations, which supports our analysis, and might be 

Ch ' applicable to other systems. 
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1 Introduction 

The dynamics of a heteropolyiner chain is relevant for the problem of pro- 
tein folding, and also from a fundamental point of view. Since a protein 
is composed of monomers of different chemical natures, it is important to 
understand the effect of heterogeneity on the kinetics of a polymer chain. 
Such results might give insight into the possible folding pathways of proteins 
or longer chains. Although the influence of quenched or annealed disorder 
on the thermodynamics of polymers is a largely investigated area of study, 
[H] together with the effect of random fields 0, random charges along the 
chain (polyelectrolytes and polyampholytes) [^, in solution or at the interface 
between two fluids |^, little is known about dynamics. 

Previous studies concerning the statics of heteropolymers show that there 
exists a frozen phase at low temperature, very similar to a spin-glass phase 
[|], which is a non-ergodic state characterized by a very slow relaxation. 
Recent phenomenological and analytical developments have reproduced the 
experimental evidence of aging in spin-glasses 0. Below a certain temper- 
ature, the system relaxes in a slower and slower way as the waiting time 
-which is the time elapsed between the beginning of the experiment and the 
observation time- is increased; the dependance on the waiting time is clearly 
proved experimentally in spin-glasses and in other glasses 0, pj. The re- 
laxation follows a power law and both time-translational invariance and the 
fluctuation dissipation theorem are violated. Similar properties have been 
found theoretically for the study of large-time out of equilibrium dynamics 
of a manifold in a random potential 0. These results are of great interest 
for systems where disorder or frustration are present, and similar ideas start 
to be applied, for instance for the rheology of soft glassy materials [0, or 
for the dynamics of structural glasses, where experimental evidence of the 
violation of the fluctuation-dissipation theorem has recently been reported 

Concerning the dynamics of heteropolymers, few studies exist at present 
[O and some of them show that there may be some glassy behaviour as the 



temperature is lowered [T^, or that the relaxation should follow a stretched 
exponential law [jl4[. There is indeed numerical evidence of stretched expo- 
nential relaxation for randomly branched polymers [ll^], or for the reptation 



of polymers in disordered media [jT6[, as well as experimental evidence of the 
same phenomenon for glasses and proteins [0 . Moreover, there is a growing 



literature about the dynamics of homogeneous but strongly frustrated poly- 
mer systems, such as polymer melts where stretched-exponential laws are 
observed through computer simulations [|T^ , and with dynamics very similar 



to the one observed in structural glasses or supercooled liquids [jT9[. Other 
approaches to the problem of dynamics of heteropolymers in melts use the 
concept of reptation to compute in a phenomenological way the relaxation 
time of a disordered chain pO |. 



Our work concerns the study of the Langevin dynamics of a heteropoly- 
mer, treated in the Mode-Coupling Approximation scheme (MCA). Our mo- 
tivation for using this approximation is based on several previous studies that 
led to significant results. This procedure goes beyond perturbation theory 
(though not in a very controlled way) and is therefore useful when one wants 
to study strong-coupling effects. The scheme is to expand the microscopic 
quantities involved in the Langevin equation to lowest non-trivial order in 
the potential -as if it was a perturbation procedure-, and then to replace in 
the correction terms the bare correlation functions (those corresponding to 
the problem without potential) by the full correlation functions one wants 
to compute. This amounts to resumming a certain class of diagrams and 
hence to go further than the weak-coupling regime. This type of procedure 
has been used by Kraichnan in the context of turbulence as a way to find 
the Kolmogorov laws starting from the Navier-Stokes equation |^; it has 



also been used for the KPZ equation where exponents close to those found by 
dynamic renormalization group were computed ||2^. Interestingly, it has also 



been found that the MCA is exact for some special models with quenched 
disorder which dynamics can be studied exactly in a mean-field approach, 
using functional methods; this is true in particular for the p-spin spherical 
spin-glass model p4|. So one can hope that the MCA approach is able to 



capture dynamic effects such as aging that arise from the presence of disorder 
in a non-perturbative way. Finally it has been pointed out that the general 
coupled set of equations obtained through MCA looks very much like those 
found in the context of the mode-coupling theory introduced by Gotze [|1^ 



which gives a reliable description of the slow dynamics of supercooled liquids, 
reinforcing the link between glassy systems which are frustrated but contain 
no disorder, and disordered systems such as spin-glasses. 

In the following, we show that the same approximation (MCA) can be 
used for the dynamics of a disordered polymer and that out-of-equilibrium 
features can be found as well. These results may be of relevance for het- 



eropolyiner melts, or for very long chains of heteropolymers. From a protein- 
oriented point of view, such results may not be directly applicable, since 
they are derived for an infinite and purely random system. However, for 
large proteins, one may observe some intermediate slow regime of folding in 
the globular state, between the fast initial hydrophobic collapse, and the final 
relaxation towards the native state once a nucleus pSl has been formed. 



2 The mode-coupling approach 

2.1 Formal developments 

We introduce here a model of heteropolymer dynamics and explain how to 
derive a set of coupled integro-differential equations using the Mode Coupling 
Approximation. We use a standard hamiltonian for a disordered polymer, 
where a quenched potential V{s, (f){s,t)) is applied and comes from the ran- 
dom nature of the interactions between monomers. In our notations, (p{s,t) 
is the position of monomer s at time t, s being the coordinate of the monomer 
along the chain, s = 1, .., A^. d is the dimensionality of the space and oq is 
the Kuhn length. 

n^±,Jis±{?^)\jisV(sMs,t)) (1) 

More explicitly, the random potential is: 



V{s, (f){s, t)) = / ds'B{s, s')6{(j){s, t) - (j){s', t)) 

We use a bar to perform the average over the disorder, and we assume that 
the value of the interaction between too different monomers is distributed in 
a gaussian way. 



B[s,s') = 



B(si, s;)5(s2, 4) = Bl5{si - S2)5{s[ - 4: 



We consider the Langevin equation for such a polymer: 
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dt Gq ds"^ 



<Pais,t) 



d 



')a{s,t) 



dsy(s,0(s,t)) + ?7,(5,t) (3) 



with a gaussian thermal noise r/(s, t), 



<Va{s,t) >=0 

< r]a{s, t)rip{s\ t') >= 2T6{s - s')5{t - t')6ap 

Our aim is to compute correlation and response functions, or at least know 
their qualitative behaviour with time. Following earlier studies, we don't 
assume a priori time-tranlational invariance and we define respectively the 
correlation function and the response function as quantities depending on 
two distinct times, t and t' . 






(4) 



Gis,t;s',t') = -Y: 
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a=\ 



d < (l)a{s,t) > 



\^J2<Ms,t)Vais',t')> (5) 



The last identity holds as long as the random noise ff{s,t) is gaussian. 
During all this study, we shall use Fourier transforms, which we define, both 
for the position 0(s, t) and for the correlation functions, with u>n = ^- 



^"(t,^ 



N Jo 



N 



Ju)n(s'-s) 



G{s, t; s',t')dsds' 



The standard procedure in the MCA is to find the solution (f)^{t) of the 
Langevin equation, up to the first non-zero order in perturbation. The dy- 
namic equation can be rewritten in Fourier space: 



dt 



<€{t)-xw:{t)+ra{t) 



(6) 



where we have added for convenience the coefficient A as the perturbative 
parameter; A is eventually set back to 1 at the end of the expansion. And 
the quantity W^{t) is defined as: 



W^it) 



N Jo 



N 



:,l^nS. 



N 



ds / ds'Bis, s') 



dq 



-iq g«Z]c«'?«('^«('''*)~'^«(*''*)) (7) 



If A is equal to we are reduced to the 'bare' problem of an elastic chain in a 
harmonic potential. Then the Langevin equation is exactly solvable and the 
solution is 



CoW 



dt'G-{t,t')f,-{t) 



(8) 



where GQ{t,t') is the bare response function, and: 



,-^l(t-t') 



G'^(t,t') = e-""^*-*'^ (9) 

When one now adds the disorder-dependant term W^{t) in the Langevin 
equation, then 



€it) 



dt'G-it,t') f,-(t') - xw:it') 



(10) 



is the exact solution, which actually gives an implicit equation for the quan- 
tity (pa{t) that cannot be solved in a straightforward way. One then performs 
in this expression an expansion up to second order in A, we refer the reader 
to Appendix A for more details. 

It is then rather straightforward, though computationally lenghty, to com- 
pute G"(t,t') and C^{t,t') as functions of the bare quantities, which are at 
the end replaced by the full or 'renormalized' quantities. One finally ends 
up with a set of coupled equations which solutions will in principle fully de- 
scribe the dynamics of the original system. We can write these equations in 
a compact way: 



d ul 
dt Oq 



G«(t, t') = 5{t -t')+ / dtiRnit, ti)&\t, t') 



+ S„(t,ti)G"(ti,t') 
Jt' 



(11) 



d u\ 

dt Qq 

t 



Jo Jo 



+ rftiS„(t,ti)C"(ti,t') (12) 

Jo 

All quantities Rn{ti,t2), S„(ti,t2), Dn(ti,t2) involved in the coupled set 
of equations are defined in Appendix A, and depend only on G'^(t, t') and 
CP{t,t'), with p ^ 0. Similar sets of coupled equations have been already 



encounted for example in [26|, p^, M and can either be solved numerically 



or require additional assumptions to get more information on the solutions. 

2.2 Analysis of the equations 

A very difficult task is to solve the set of integro-differential equations de- 
scribed above. One major difficulty relies in the fact that all modes are 
coupled, and as it has already been observed for example for the mode- 



coupling equations of the KPZ model [^ |2^, the numerical treatment for 
these equations presents numerous problems. We have not been able to make 
significant progress in that direction; not only should it require a recursive 
algorithm with careful check for the convergence of all functions, but we also 
expect some divergences in the long time regime, which would require the 
introduction of unknown cut-offs. 

In our analysis, we took into account all terms found in the former section, 
without truncating them with too crude approximations. However it would 
be interesting in the future to find a way to simplify these equations (even if 
the connection with the initial model becomes then less obvious), that would 
reproduce the results that we find here. 

The first step in the analysis can be done by focusing on the large time limit 
and proposing an ansatz for the correlation functions in that time domain. 
Let's assume that one can write, in the limit where t — *> oo and t' -^ oo, with 



C''it,t') = qUl-\ (13) 

mt,t') = ^fJ^y (14) 



Such an ansatz also contains a generalized version of the fluctuation-dissipation 
theorem (this has already been introduced in earlier studies [|^), which can 
be written: 

G^{t,t')=x^C''{t,t') (15) 

where x is the coefficient x = —. 

■yq 

Then by replacing into the mode-coupling equations, one is left with implicit 
equations for q, f^, x and 7, but the dependance on time cancels out in the 
limit of large times. This is what makes the ansatz consistent, at least as 
far as the dependance on time is concerned. We give computational details 
in Appendix B. In particular, in the equation for C"(t,t'), D„(ti,t2) can be 
written as a sum of four contributions D^\ti,t2), z = 0, 1,2,3 as shown in 
Appendix A, and the terms involving D^\ti,t2) and D'^\ti^t2) are shown 
to be negligible in the limit of large times; this makes the time cancellation 
possible. All parameters g, /„, x and 7 have a dependance on temperature. 
Then the remaining equations on q, fn, x and 7 are again too difficult to solve 
numerically, since they require the introduction of cut-offs - see Appendix B 
-, to stop the divergences in the integrals when the two time arguments 
become too close to each other (in particular, both G^{t,t') and C^{t,t') 
have different analytical forms when t' —>■ t, which we don't know). 

The two extreme cases T = and T = 00 are of interest in this problem. 
In the case where the temperature is zero, the equations can be simplified, 
the terms involving D^\ti, ^2) and D^^)(ti, ^2) are equal to zero and one sees 
easily that the above ansatz still remains valid, for the same reasons as the 
ones explained above (see Appendix B). 

If we study the limit T -^ 00, one is left with a single term in the equation 
for ^"(t,^): 



/* t/tiDf (t,ti)G'"(t',ti) = (16) 

Jo 



The power-law ansatz used above can no longer satisfy this condition. 
However, if one assumes time-tranlational invariance and the usual fluctuation- 
dissipation theorem (G(r) = — y'^tC{t)) the mode-coupling equations are 
simply satisfied. In the limit of infinite temperature, one is actually left with 
the simple Rouse model for a homopolymer chain, and an exponential time 
relaxation. 



Such information about a glassy behaviour at low temperatures with 
power-law behaviour lead naturally to think of the existence of a critical 
temperature T^ that may separate the glassy, non-FDT, non-TTI regime 
from a high temperature regime where the relaxation would be typically ex- 
ponential. Although it is not possible to determine this temperature from 
our equations, it should be easier to observe such a phenomenon in simpler 
models of polymers, or by studying numerical models of polymers. 

3 An alternative treatment of the mode-coupling 
equations: Functional renormalization group 
approach 

In view of the difficulties raised by the mode-coupling equations, one has to 
search for new analytical methods to try and solve them. One of them is 
to apply a functional renormalization group analysis to the mode-coupling 
equations themselves. To our knowledge, such a method has never been used 
in this context. For the present problem, this procedure enables us to have 
more information about the analytical form of the correlation functions. In 
particular it can justify some scaling form for their analytical expressions, 
as soon as one finds a fixed point in the RG procedure that is believed to 
represent the small frequency, small wavevector regime. This is motivated 
by the fact that we are mostly interested in the long times limit, and in the 
long-distance regime {s — s' ^ oo, along the chain). 

In the case of the disordered polymer, our dynamic RG calculation gives 
rise to a fixed point. However, the fixed points equations are themselves hard 
to solve. We also believe that such a method could be of interest for simpler 
and largely studied systems, such as the KPZ equation. 



In this calculation ||2^, we want to integrate out the fast wavevector 
modes and keep only the slow modes, the ones with small wavevectors. One 
also has to do the same for high frequencies, to keep track only of the low 
frequency part; this can be done by expressing all the quantities in frequency 
space, but we didn't describe it here for simplicity. For that purpose, we 
first switch from discrete to continuous Fourier variables, by replacing J2n by 
/ ^ with an upper cut-off A. The wave vectors such that ^ < ||A;|| < A can 
be integrated out, with b = e^^ being close to 1; then the only perturbative 



parameter here is 51. After integration, we denote the quantities Q for which 
the ||A;|| integration is only /q ' ^ by Q^. 

In a more convenient way than the ones used in the previous section, the 



starting mode-couphng equations we used can be rewritten in the form p4 
(see also Appendix A): 



Ck{t,t') = Cl{t,t') + f dh f dhGk{tM)Dk{tiM)Gk{t\h) (17) 

JO Jo 



Gk{t,t') = Gl{t,t')+ f dhGUt,h)\ r dt2Rk{ti,t2)Gk{h,t') + r dt2T.,{t,M)Gk{t2.t') 
Jo Uo Jo 

(18) 

In such a way, we can write for a quantity Qk{ti,t2), the following one- 
order expansion {Q = G,C, R,T,, D): 

Qk{ti,t2) = Qt{ti,t2) - 6lQk{h,h) (19) 

The corrections Qk{ti,t2) can be computed and only involve the different 
quantities Gq and Q^. In particular, for Qkit^t') and Ck{t,t'), we have: 

ft i-ti rti 



Gk{t,t')= f dhGl{t,h)[f ' dt2'JZk{h,t2)Gk{ti,t') + r dt2Sk{ti,t2)Gk{t2..t') 

Jo Jo Jo 

+ / ' dt2Rk{tl,t2)gk{tl, t')+ r dt2'Lk{h,t2)Qk{t2, t')] (20) 

Jo Jo 



Ckit,t')= f dh I dt2[Gk{t,ti)Dk{tiM)Gk{t'M) 
Jo Jo 

+Gk{t, h)Vk{h, t2)Gk{t', t2) + Gk{t, h)Dkih,t2)Gkit', ta)] (21) 

The next step is to do a rescaling of all quantities and write a differential 
equation for the renormalization flow where the increment is 61. One has 
then to make scaling assumptions, which are expected, at least at the fixed 
points, if any. Then, if one assumes 

Qkit,t') = ^qitk\t'k^) (22) 

10 



where z is the dynamic exponent, the renormahzed quantity is 

Qf(t,t') = 6>^Q,l(r%r^0 (23) 

By expanding this last expression to first order in bl one finally ends up with: 

More specifically we used for the q functions the scaling forms that correspond 
to the ansatz of the preceding section (see equations (13) and (14)), written 
in such a way that the values of x can be simply identified, 

C.("') = p^': (j) (25) 

'''<*•''' -F^IF'O '''^ 

O'.C-O = p^d (0 (27) 

Due to the scaling nature of all quantities, the derivative terms -S^ and 
-^ in the fiow equations can be simplified and expressed in terms of Qk{t, t'), 
and we end up with the following set: 



dl 

dCk 
dl 

dPk 
dl 



zGkit,t')+gkit,t') (30) 

-2{d + 2a-z)Ck{t,t')+Ck{t,t') (31) 

-2{d + 2a- 3z)Dkit, f) + V^it, t') (32) 
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-^ = zRk{t,t')+nk{t,t') (33) 

^ = zJ:k{t,t')+Sk{t,t') (34) 

The fixed points are obtained by setting the / derivative to 0, and denoting 
by Q* the fixed-point quantities, and replacing in the expressions for Qk{t, t') 
and Ck{t,t'), we obtain the self-consistent equations at the fixed point. 



Gl{t, t') = f dt.Glit, t^) \ r dt2Rl{h, h)Gl{h, t') + r dt2El{t,,t2)Glit2, t') 
Jo uo Jo 

(35) 

{d + 2a-z)Cl{t,t') = {d + 2a-Az) [ dh [ dt2Gl{t,ti)Dl{ti,t2)Gl{t' ,12) 

Jo Jo 

(36) 

We checked that the power-law ansatz used in the first section is still 

a solution for these equations of fixed point; the procedure of replacing the 

ansatz in the expressions is exactly the same as the one described in Appendix 

B. This can then justify its use in the mode-coupling equations (11) and 

(12) from which we started our analysis, and which represent the real-space 

quantites of interest; the scaling form of the ansatz is also justified. If we 

then plug at the fixed point the generalized fiuct nation- dissipation theorem 

G^{t,t') = x^C^{t,t'), we obtain the following relation between exponents: 

2a+d = 2z. Further information on the exponents a and z could be obtained 

by solving numerically coupled equations like (35) and (36); this will be done 

elsewhere. 



4 Discussion and Conclusion 

The analysis of the mode-coupling equations using renormalization group 
ideas enabled us to write the dynamic correlation functions of a disordered 
polymer in terms of scaling functions. Unfortunately, at this point, it does not 
give much insight into the exact asymptotic time and wavevector dependance, 
at large times and large wavelenghts. We believe however that this new 
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method could be very convenient for more simple examples of stochastic 
equations, that are treated using mode-coupling techniques. 
The main result of the mode-coupling approach is the evidence of the out- 
of-equilibrium character of the dynamics of heteropolymers, in the thermo- 
dynamic limit. The dependence of the dynamic correlation functions on two 
times -the smaller one being the waiting time-, implying no time-translational 
invariance and a generalized fluctuation-dissipation theorem, is required at 
low temperatures; more precisely, a power-law dependence on two times is a 
solution of our set of equations. This contains the aging phenomenon, similar 
to the one observed in spin-glasses and other types of disordered systems. 
Similar calculations as the one described above can be made, by slightly 
changing the starting Langevin equation. We replaced the disordered in- 
teraction B{s,s') between monomers by a constant attractive interaction 
coefficient v, therefore restricting our study to homopolymers. Remarkably 
enough, we found that the power-law ansatz still holds at low temperature. 
Even in the absence of disorder, the mode-coupling equations seem to induce 
some apparent disorder in frustrated systems (such as polymer chains in the 
collapsed phase in this case), and lead to non-trivial glassy behaviour. In the 
MCA approach, it appears that one can capture the frustration character of 
a collapsed polymeric chain that stems from the competition between the at- 
tractive interaction between monomers and the harmonic potential due to the 
elasticity of the chain (the classical Mode Coupling Theory [|1^] for glasses 



also captures frustration in such simple realizations as binary mixtures of 
Lennard- Jones -LJ- particles). 

This similarity of behaviour between disordered and non-disordered polymers 
is not very surprising when one remembers that mode-coupling equations 
have very similar structures in the disordered and non-disordered version 



of a given model [Q. More precisely, it has been previously shown that 
given a non-disordered but frustrated or chaotic model, one can find a disor- 
dered model of the same class whose mode-coupling equations are actually 
exact, and coincide with the approximate mode-coupling equations of the 
non-disordered model. (This is however not exactly the situation that occurs 
here, since the mode-coupling closure is only an approximation scheme for 
our heteropolymer model.) This emphasizes the importance of such studies 
for real glasses, which are not characterized by the presence of disorder. In 
particular, our analysis may be helpful to understand the dynamics of poly- 
meric glasses, or melts, that also exhibit aging [Q] and vitreous dynamics 
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1^, qualitatively very similar to the physics of binary Lennard- Jones fluid 



mixtures 28 



It is interesting to note that a great similarity has been observed in the dy- 
namics of supercooled liquids ( modelled by binary mixtures LJ particles) 
and of polymer melts, whereas a few differences are noticeable due to the 
connectivity effect of the polymeric system [|18[. It has been shown that 
the correlation functions and in particular the dynamic structure factor for 
polymer melts can be well described by the MCT for liquids. However the 
mean-square deplacement shows a subdiffusive behaviour compared to the 
case of non-connected LJ particles, due to the connectivity of the chains. 
This difference occurs in the a regime, at large times, that reflects the rear- 
rangements of the chain at scales larger that the size of a particule's cage, 
whereas the /3 regime, at small times, corresponds to local movements inside 
a cage, which are the same in the presence or not of connectivity. On the 
ground of these microscopic mechanisms, it does not seem surprising either, 
that a disordered polymer melt or a homogeneous melt may have the same 
kind of dynamics. 

To summarize, our model provides a tentative theoretical framework which 
allows to exhibit slow dynamics and aging in dense phases of homogeneous or 
disordered polymers. A further step in this study would be to find an efficient 
numerical algorithm for the resolution of the mode-coupling equations or of 
some simplified form of these mode-coupling equations. This would provide 
the full solutions for the correlation functions, as well as for the mean-square 
displacement, as functions of time, and allow for a direct comparison with 



the numerical results of Baschnagel and al. [|T8|. 

An important issue that also needs to be clarified is in which way the 
fiuctuation-dissipation theorem is violated at long times. It has been shown 
very recently |^ that in the case of systems with short range interactions, the 
measurement of the violation of the dissipation-fiuctuation ratio is directly 
related to the pattern of the replica symmetry breaking ('one-step rsb' or 'full 
rsb'). The statics of the heteropolymer model studied here were reported in 
[H] and showed a one-step replica symmetry breaking, so the fiuctuation- 



dissipation should, according to ||29| be of the form given by equation (15). 



In addition, concerning the case of homopolymers, it would be of great in- 



terest to use the first principle computations described in |^ and be able 
to draw a link between the equilibrium and the dynamical properties in the 
glass phase. 
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Appendix A 

In this appendix, we give the full expressions for the mode-coupling equa- 
tions obtained in section 2, as well as some intermediate results. 

Starting from equation (10), one performs an expansion to second order 
in A, which leads to an approximate expression for 0"(t). This expression is 
the following: 



2A r* 

+ —- / dt'Gm.t') / ds / ds'e'^-'B{s,s') 
N Jo Jo Jo 



Y. ^qp j:i^-'-'' - <^~'^'"') Jl dt"Gl{t\ t")W$^,{t") 



fS p^O 



where 



Starting from this expression we calculated G"(t, t') and C"(t, t'), performing 
both averages on the random noise and on the disorder. The result can be 
written in terms of 0^(1,1') and C^it^t') only; but the MCA approximation 
consists in replacing these quantities by the full unknown quantities G'^{t, t') 
and CP{t, f), resulting in a set of coupled Dyson equations: 



Cnit^t') = Cl{tX) + f dh f dt2Gn(t,t^)D^{hM)Gn{t' M) 

Jo Jo 



Gn{tX) = Gl{tX)+ f dhGl{tM) \r dt2Rn{hM)Gn{hX) + T rft2S„(ti, t2)G„(t2, t') 

JO iJo Jo 

The mode-coupling equations can be rewritten as integro-differential equa- 
tions, performing on the left and right-hand sides of the previous equations a 
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^ -1 



convolution with Go • The resulting set of equations is the one introduced 
in section 2: 



dt Oq 



G'\t, t') = d{t -t')+ f dtiRnit, ti)G"(t, t' 

Jo 



+ S„(t,ti)G"(ti,t') 
Jt' 



d ujI 
dt Oq 



C''it,t') 



*(itii?„(t,ti)c'"(t,t')+ r rftiz}„(t,ti)G'"(t',ti) 



+ / rftiS„(t,ti)C'"(ti,0 



The term /q dtiRnit, ti) plays the role of a mass and is given by the following 
expression: 



Rn{tl,t2 






f-N f-N 

/ dsi / ds[{l - e'-"(4-^i)) y (1 - e"^^^'''<^)GP{ti,t2 
Jo Jo ^0 



[A(si,s;;ti,t2)] ' -^3Xis^,s[■,t,)Yis^,s[■,h,t2) 

The quantity S„(i(:i,t2) can be seen as the self-energy associated with the 
corresponding Dyson equations: 



Sra(ti, ^2 



2S 



A^3 (^2nY Jo Jo 



dsi I ds\{l - e^'^"("i-^i)) y (1 - e^^^("i'<))GP(ti,t2 



[A(S1, S\- ti, t2)] ' ^ [X(S1, S'l; ti)X(si, S[- t2) + 2F2(s,, g'^; ti, t2) 

Finally, D„(ti,t2) can be seen as a 'renormalized' noise correlator, which 
expression is reported here: 

Each term can be computed separately: 
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IT 



ABl 1 fN 



[A(si,s;;ti,t2)] 



/ dsi / (is'^y(si, s'^; ti, ^2 
JO Jo 



i-i-f 



„/'^^ / N 2i3n 2T 1 






iV2 N {27iY Jo ' Jo 
^ /'' dt.G^it^, t3)(l - e*-"(^i-<)) [A(si, s'l; ti, t2 



PT^O 



;-3X(si, s'l; t3)F(si, s[; t^, h) / rft4(5(ti - t4)G"(t2, ^4 



t2 



+ 



X(si, s'l; t2)^(si, 4; ^3) + 2F2(si, 4; ^2, ts) / dU5{ti - t4)G'"(t3, ^4)) 



Finally we used the following notations: 

017^0 

niy^O 



A{si,s[;ti,t2) = X {si, s[; ti)X {si, s[; t2) -Y'^{si,s[;ti,t2) 



Appendix B 

In this appendix, we show how the ansatz proposed by equations (13) 
and (14) can be a suitable solution for the mode-coupling equations obtained 
in section 2. 

We give the expressions of the different terms introduced in Appendix A, 
after performing adequate change of variables, after the power-law ansatz has 
been pluged into these terms. Whenever required, we introduced lower and 
upper cut-offs, assuming that in the extreme regimes of times (short times 
or t' -^ t, for which we don't have explicit analytical forms) the integrals 
are actually convergent, and that the time scaling we find is not modified by 
these contributions. In all equations we neglect the derivative terms, which 
is justified in the limit f <^ t. 

!• Equation for Gn{t, t') 

• (i) The first term in front of the propagator in the right-hand side of 
the equation is comparable to a mass: 

t IB"^ ^ r^ 3(7' 



.„ *ft.(MO = -^p^y„ "s(l - e— )^|A-(. 



'"d.m-'+^-'d -«=')-'-« 



This factor is independant of time, up to the cut-off 1 — e. 
(ii) Concerning the second term, it is in fact proportional to G'^{t,t'): 
rt . 27?2 1 /-AT a''^f 



]" duu-\l - u^--)-^-^i[l + u^-<\j {^-\ 



Finally, given the cut-offs w and 1 — e, all terms in the equation for 
G'"(t, t') are proportional to each other, and hence the dependance on 
time cancels out. 

2« Equation for C'„(t,t') 
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(i) The first new term that has to be computed is in fact proportional 
toC'"(t,t'): 

t ^ 9R2 1 /-AT „/f 



n _ zH r^^ n t 



(ii) The remaining terms coming from the -D„(t, ti) contribution can be 
simphfied. First, when t' < t, it is easy to see that the ones involving 
-D^°)(t,ti) and D^^\t,ti) are equally zero. 

(iii) The D^\t,ti) contribution is again proportional to C'^{t,t'), 



1-d 



(iv) Finally, the last term has a different time dependance. 

ft' 9 R2t^ 1 ^7V „/ 

f-3 f'duu"'-^ f''^dvv^^-'[l~v'T'~~' + f'duu"'-^ f'~^dvv-'[l-v'T''^A 

\ Jw Jo Jw Jw / 

Because of its scaling behaviour in t and t', this term is in fact a sub- 
dominant term compared to C"'{t, t') in the limit t — > cxd and f — >■ oo, 
t' <^ t, and we can neglect it. 



20 



References 

[1] T. Garel, H. Orland and E. Pitard, "Protein folding and heteropoly- 
mers", in Spin Glasses and Random Fields, A. P. Young (ed.), World 
Scientific, Singapore (1997). 



[2] 
[3] 

[4] 
[5] 

[6: 



[7; 

[9 
[10 

[11 
[12: 



S.F. Edwards, M. Mutliukumar, J. Chem. Phys., 89, 2435 (1988). 

M. Kardar, Y. Kantov, "Randomly charged polymers" in Complex Be- 
haviour of Glassy Systems, M. Rubi, C. Perez- Vicente (eds.). Lecture 
Notes in Physics 492, Springer- Verlag (1997). 

T. Garel, D.A. Huse, L. Leibler, H. Orland, Europhys. Lett, 8, 9 (1989). 
C Monthus, preprint |cond-mat/9902336 . 



T. Garel, H. Orland, Europhys. Lett., 6, 307 (1988). E.I. Shakhnovich, 
A.M. Gutin, J. Phys. A, 22, 1647 (1989), Europhys. Lett, 8, 327 (1989). 

J-Ph. Bouchaud, L.F. Cugliandolo, J. Kurchan, M. Mezard, Out of equi- 
librium dynamics in Spin Glasses and other Glassy Systems in Spin 
Glasses and Random Fields, A. P. Young (ed.). World Scientific, Singa- 
pore (1997). 

M. Alba, M. Ocio, M. Hammann, Europhys. Lett, 2, 45 (1986). 

L.C.E Struik, Physical Aging in Amorphous Polymers and Other Mate- 
rials, (Elsevier, Houston, 1978). 

L.F. Cughandolo, P. Le Doussal, Phys. Rev. E, 53, 1525 (1996), and 
references therein. L.F. Cugliandolo, J. Kurchan, P. Le Doussal, Phys. 
Rev. Lett, 76, 2390 (1996). 

P. Solhch, F. Lequeux, P. Hebraud, M.E. Gates, Phys. Rev. Lett., 78, 
2020 (1997). P. Sollich, preprint |cond-mat/97l200l| . P. Hebraud, F. 



Lequeux, preprint |cond-mat/ 9805373 . 



T.S. Grigera, N.E. Israeloff, preprint |cond-mat/990435l| 



J-R. Roan, E. Shakhnovich, Phys. Rev. E, 54, 5340 (1996). 



21 



[13] D. Thirumalai, V. Ashwin, J.K. Bhattacharjee, Phys. Rev. Lett., 77, 
5385 (1996). 

[14] M. Rehkopf, V.G. Rostiashvili, T.A. Vilgis, J. Phys. II France, 7, 1469 
(1997). V.G. Rostiashvili, M. Rehkopf, T.A. Vilgis, Eur. Phys. J. B, 6, 
233, 497 (1998). 

[15] J.R Kemp, Z.Y. Chen, Phys. Rev. E, 56, 7017 (1997), Phys. Rev. E, 
60, 2994 (1999). 

[16] D. Cule, T. Hwa, Phys. Rev. Lett, 80, 3145 (1998). 

[17] S.J. Hagen, W.A. Eaton, J. Chem. Phys., 104, 3395 (1996). J.L. Green, 
J. Fan, C.A. Angell, J. Phys. Chem, 98, 13780 (1994). I.E.T. Iben and 
al., Phys. Rev. Lett, 62, 1916 (1989). 

[18] K. Binder, J. Baschnagel, C. Bennemann, W. Paul, J. Phys: 
Cond. Matter, 11, A47 (1999). J. Baschnagel, C. Bennemann, W. 
Paul, K. Binder, Dynamics of a supercooled polymer melt above 
the mode-coupling critical temperature: cage versus polymer- specific 
effects, and references therein, contribution to the workshop Uni- 
fying Concepts in Class Physics, Trieste, Sept 15-18, 1999, see 
|http://www.ictp.trieste.it^franz/workshop.html| , to be published in J. 
Phys: Cond. Matter. 

[19] W. Gotze, Les Houches 1990: Liquids, Freezing and the Class Transition 
part 1, ed J-P. Hansen, D. Levesque, J. Zinn- Justin (North Holland), 
287-503. W. Gotze, L. Sjogren, Rep. Prog. Phys., 55, 241 (1992). For a 
recent review, W. Gotze, J. Phys: Cond. Matter, 11, Al (1999). 

[20] J-P. Bouchaud, M.E. Gates, J. Phys. II, 3, 1171 (1993). A.N. Semenov, 
Phys. Rev. E, 60, 3076 (1999). 

[21] R. Kraichnan, J. Flmd. Mech, 5, 497 (1959). 

[22] T. Hwa, E. Frey Phys. Rev. A, 44, R7873 (1991). E. Frey, U.G. Tauber, 
T. Hwa Phys. Rev. E, 53, 4424 (1996). J.P Bouchaud, M. Gates, Phys. 
Rev. E, 47, R1455 (1993), Erratum 48, 635 (1993). M. Moore and al. 
Phys. Rev. Lett, 74, 4257 (1995). 



22 



[23] Y. Tu, Phys. Rev. Lett 73 3109 (1994). 

[24] J.P Bouchaud, L. Cugliandolo, J. Kurchan, M. Mezard, Physica A, 226, 
243 (1996). 

[25] V.I. Abkevich, A.M. Gutin, E.I. Shakhnovich, Biochemistry, 33, 10026 
(1994). 

[26] L.F. Cugliandolo, J. Kurchan, Phys. Rev. Lett, 71, 173 (1993). 

[27] P.C. Hohenberg, B.I. Halperin, Rev. Mod. Phys. 49 435 (1977). 

[28] W. Kob, J-L. Barrat, Physica A 263 234 (1999) and references therein. 

[29] S. Franz, M. Mezard, G. Parisi, L. Peliti, Phys. Rev. Lett, 81, 1758 
(1998); preprint |cond-mat/990337q . 



[30] M. Mezard, G. Parisi, A first principle computation of the thermody- 
namics of glasses, and references therein, contribution to the work- 
shop Unifying Concepts in Glass Physics, Trieste, Sept 15-18, 1999, see 
|http://www.ictp.trieste.it^franz/workshop.html| , to be published in J. 
Phys: Cond. Matter. 



23 



